SYMPLECTIC SEMICLASSICAL WAVE PACKET DYNAMICS 
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Abstract. The paper gives a symplectic-geometric account of semiclassical Gaussian wave packet 
dynamics. We employ geometric techniques to "strip away" the symplectic structure behind the 
time-dependent Schrodinger equation and incorporate it into semiclassical wave packet dynamics. 
We show that the Gaussian wave packet dynamics of Heller is a Hamiltonian system with respect 
to the symplectic structure, apply the theory of symplectic reduction and reconstruction to the 
dynamics, and discuss dynamic and geometric phases in semiclassical mechanics. We also propose 
an asymptotic approximation of the potential term that provides a practical semiclassical correction 
term to the approximation by Heller. A simple harmonic oscillator example is worked out to 
illustrate the results, along with the canonical and action-angle coordinates for the system. Finally, 
we look into the geometry behind the Hagedorn parametrization of Gaussian wave packet dynamics. 



1. Introduction 

1.1. Background. Gaussian wave packet dynamics is an essential example in time-dependent 
semiclassical mechanics that nicely illustrates the classical-quantum correspondence, as well as a 
widely-used tool in simulations of semiclassical mechanics, particularly in chemical physics (see, 
e.g., Tannor [39j and Lubich [24J). A Gaussian wave packet is a particular form of wave function 
whose motion is governed by a trajectory of a classical "particle"; hence it provides an explicit 
connection between classical and quantum dynamics by placing "(quantum mechanical) wave flesh 
on classical bones." [6^ 39j 

The most remarkable feature of Gaussian wave packet dynamics is that, for quadratic potentials, 
the Gaussian wave packet is known to give an exact solution of the Schrodinger equation if and only 
if the underlying "particle" dynamics satisfies a certain set of ordinary differential equations. Even 
with non-quadratic potentials, Gaussian wave packet dynamics is an effective tool to approximate 
the full quantum dynamics, as demonstrated by, among others, a series of works by Heller |16| 
\T7\ [T8] . See also Russo and Smereka [35] for a use of the Gaussian wave packets to transform the 
Schrodinger equation into more computationally tractable equations in the semiclassical regime. 

It also turns out that Gaussian wave packet dynamics has nice geometric structures associated 
with it. Anandan showed that the frozen Gaussian wave packet dynamics inherits symplectic 

and Riemannian structures from quantum mechanics. Faou and Lubich |10j (see also Lubich [2^ 
Section II. 4]) found the symplectic/Poisson structure of the "thawed" spherical Gaussian wave 
packet dynamics (which is more general than the frozen one) and developed a numerical integrator 
that preserve the geometric structure. It is worth noting that Heller [TBJ decouples the classical 
and quantum parts of the dynamics and only recognizes the classical part as a Hamiltonian system, 
whereas Faou and Lubich [10] show that the whole system is Hamiltonian. 

1.2. Main Results and Outline. The main contribution of the present paper is to provide a 
symplectic and Hamiltonian view of Gaussian wave packet dynamics. Our main source of inspiration 
is the series of works by Lubich and his collaborators compiled in Lubich |24j. Much of the work 
here builds on or gives an alternative view of their results. Our focus here is the symplectic 
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point of view, as opposed to the mainly variational and Poisson ones of Faou and Lubich [TO] and 
Lubich [2i]. In Section [2| we start with a review of some key results in [24J from the symplectic 
point of view, and then consider the non-spherical Gaussian wave packet dynamics in Section |3j 
The main result in Section [3] shows that Heller's non-spherical Gaussian wave packet dynamics is 
a Hamiltonian system with respect to the symplectic structure found by a technique outlined in 
Section [2] the result is shown to specialize to the spherical case of Faou and Lubich [lOj in Section [sj 
Then, in Section [4j we exploit the symplectic point of view to discuss the symplectic reduction of 
the non-spherical Gaussian wave packet dynamics. This naturally leads to the reconstruction 
of the full dynamics and the associated dynamic and geometric phases in Section [6j Section [7] 
gives an asymptotic analysis of the potential terms present in the Hamiltonian formulation. The 
potential terms usually cannot be evaluated analytically and one may need to approximate them 
for practical applications. We propose an asymptotic approximation that provides a correction 
term to the locally quadratic approximation of Heller. We illustrate the results by considering a 
simple one-dimensional harmonic oscillator in Section [Sj The semiclassical harmonic oscillator is 
completely integrable: We find action-angle coordinates using the Darboux coordinates found in 
Section [5] and the associated Hamilton-Jacobi equation. We also find the explicit formula for the 
reconstruction phase. Finally, in Section [9], we exploit the geometry of the Siegel upper plane to 
give a geometric interpretation of the Hagedorn parametrization of Gaussian wave packet dynamics. 



2. Symplectic Model Reduction for Quantum Mechanics 

This section shows how one may reduce an infinite-dimensional quantum dynamics to a finite- 
dimensional semiclassical dynamics from the symplectic- geometric point of view. It will also be 
shown that the finite-dimensional dynamics defined below is optimal in the sense described in 
Section 2.3 We follow Lubich [2A\ Chapter II] with more emphasis on the geometric aspects to 
better understand the geometry behind the model reduction. 

2.1. Symplectic View of the Schrodinger Equation. Let T-L he a complex (often infinite- 
dimensional) Hilbert space equipped with a (right-linear) inner product (•,•). It is well-known (see, 
e.g., Marsden and Ratiu |i26^ Section 2.2]) that the two- form onTi defined by 

is a symplectic form, and hence is a symplectic vector space. One may also define the one-form 
6 on 7^ by 

G(V') = -^Im (-0, dxP) ; (G(V'), = -him {iP, ip) . 

Then, one has Q, = — d0. Now, given a Hamiltonian operatoi]^ H : Ti ^ H, we may write the 
expectation value of the Hamiltonian {H) : "H — )• M as 

Then, the corresponding Hamiltonian flow 



on Ti defined by 

gives the Schrodinger equation 



ix^.n = d{H) (1) 



h 



In general, the Hamiltonian operator H may not be defined on the whole Ji. 
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2.2. Symplectic Model Reduction. Let be a /inzte-dimensional manifold and suppose there 
exists an embedding l : Ai ^ T-L and hence is a submanifold of T-L. 

Proposition 2.1 (Lubich [241 Section II. 1]). If the manifold A4 is equipped with an almost complex 
structure Jy : TyM — )• TyAi such that 

Tyi O Jy = i ■ Tyi (2) 

for any y £ M, then M is a symplectic manifold with symplectic form •= '•*^- 

The proof of Lubich ^24j is based on the projection from V. to the tangent space r^(j^)i(A^) of the 
embedded manifold l{A4). We give a proof from a slightly different perspective using the embedding 
I. : M. ^ H more explicitly. As we shall see later, the embedding i is the key ingredient exploited 
to define geometric structures on the semiclassical side as the pull-backs of the corresponding 
structures on the quantum side. 



Proof. It is easy to show that 0_a4 is closed: = t*dQ = 0. We then need to show that 

is non-degenerate, i.e., TyM n {TyM)^ = {0}, where (• )-*- stands for the symplectic complement 

with respect to Q.m- Let Vy G TyM. n {TyM.)-^\ then Jy{vy) E TyM. and thus 

= VL^iVy, Jy{Vy)) 

= n{TyL{Vy),TyLO Jy{Vy)) 
= 2hlm {Tyi{Vy),iTyi{Vy)) 
= 2hRe {TyL{Vy),Tyi{Vy)) 
= 2h{Tyi{Vy),Tyi{Vy)) . 

Hence Tyi{vy) = and so Vy = since l is injective. Therefore, TyM n (TyA^)-*- = {0} and thus 
M is symplectic with the symplectic form D 



Now, define a Hamiltonian : 7W — )• R by the pull-back 

H ■.= i*{H) = {H)oi. 

Then, we may define a Hamiltonian system on M by 

iXh^M = dH. (3) 

Hence we "reduced" the infinite-dimensional Hamiltonian dynamics ^(^j^'^ on T-L to the finite- 
dimensional Hamiltonian dynamics Xjj on M. 

If we write the embedding l : M ^ Ti explicitly asy xiu): then one may first find a symplectic 
one-form @m on M as the pull-back of by i, i.e., 

e^, := = -him (x, dyK (4) 

Then, the symplectic form •= i*^^ is given by 

On the other hand, one can calculate the Hamiltonian : — )• M as follows: 



H{y) = {x{y),Hx{y))- 



(5) 
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2.3. Riemannian Metrics and Least Square Approximation. As shown by Lubich [24' Sec- 
tion II. 1.2], it turns out that the the finite-dimensional dynamics Xh is the least square approxima- 
tion to the original dynamics X^^^ in the sense we will describe below. Again, Lubich [24j exploits 
the projection from Ti to the tangent space r^(j^)/,(7W), but we give an alternative account using the 
metrics naturally induced on Ti and Ai. 

First recall (see, e.g., Marsden and Ratiu [26', Section 5.3] and Chruscihski and Jamiolkowski [HI 
Section 5.1.1]) that any complex Hilbert space 7i is equipped with a Riemannian metric naturally 
induced by its inner product. In our setting, we may define 

g{ipi,ip2) ■= 2hRe {ipi,ip2) 

so that it is compatible with the symplectic structure Q in the sense that 

g{iipi,ip2) = ^{■4'i,ip2) and n{ipi,iilj2) = g{'4>i,ip2)- (6) 

Then, we may induce a metric on Ai by the pull-back 

gM ■= i^*g, 

and thus we may define norms || • || and || • for tangent vectors on Ti and Ai, respectively, as 
follows: 

ll^ll := Vs{X,X), \\v\\j^ := y^gM{v,v)- 

Proposition 2.2 (Lubich [2^1 Section II. 1.2]). // the manifold Ai is equipped with an almost 
complex structure Jy ■ TyAi — > TyAi that satisfies then the the Hamiltonian vector field Xh 
on Ai defined by Q is the least square approximation among the vector fields on A4 to the vector 
field -^(^h) defined by the Schrddinger equation For any y ^ Ai let rj := i{y) G %; then, for 
any Wy G Ty Ai, 

ll^(H>(^) - Tyi{wy)f > ||A^^^(r/)f - \\XH{y)\\l„ 
where the equality holds if and only if Wy = Xniy)- 

Proof. Notice first that the inclusion map l pulls back the compatible triple — metric, symplectic 
form, and complex structure — to Ai, i.e., Eq. ^ implies, for any v,w ^ TAi, 

gMiJiv),w) = QMiv,w) and 0^(t!, J(?i;)) = g^(t>, w). 

We may then estimate the difference between A^^^ and W := Tl[w) for any w G TAi as follows: 

\\Xi^H)-Wf = 9{x^^^-W,X^f,^-w) 

= 9 (^(H) ' ^(i/) ) - 25 (^(^) ,w)+g{W,W), 

where 

g{x^f,^,w)=^l{x^f,^,^w) 

= ^l(x^fj^,TioJ{w)) 

= d{H) ■ TioJiw) 
= d{i*{H))-J{w) 
= dH ■ J{w) 

= gM{XH,w) , 
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and g{W, W) = QMiw, w). Therefore, 

W^im - = 9\^im^X,fjA - 2gM{XH,w) + gM{w,w) 



^{H)\\ - IK^WIIX) 

where the equahty holds if and only w = Xh- When w = Xh, we obtain 

||X^^^(r,) - TyiiXHiy))f = - \\XHiy)fM, 

so the inequality can be rewritten as 

\\X^f,^{r^)-Tyiiwy)f>\\X^^^{r^)-TyL{XH{y))f, 
which is to say that Tyi{XH{y)) is the best approximation to X^^^^ (rj) in the least squares sense. □ 

3. Gaussian Wave Packet Dynamics 

3.1. Gaussian Wave Packets. In particular, let H ■= L^(M"') with the standard right-linear inner 
product (•, •) and H be the Schrodinger operator: 

H = -^A + V{x), 
2m 

where A is the Laplacian in M*^. 

Let us now consider the following specific form of x called the (non- spherical) Gaussian wave 
packet (see, e.g.. Heller p!6', T7J): 



xiy;x) = exp 



^(x - q)^C{x - q) + p ■ {x - q) + {(f) + i6) 



(7) 



where C = A + iB is a d x d complex symmetric matrix with a positive-definite imaginary part, 
i.e., the matrix C is an element in the Siegel upper plane [36| defined by 

:= |C =^ + i6 G C'^^'^ \A,Be Symrf(M), i3 > o} , (8) 

where Sym^(M) is the set of d x d real symmetric matrices, and B > means that B is positive- 
definite. It is easy to see that the (real) dimension of is d{d + 1). 
One may then let Ai be the (d + l){d + 2)-dimensional manifold 

M = T*R'^ X X X M, 

and its typical element y E is written as follows: 

y ■■= {q,p,A,B,(l),6). 

We then define an embedding oi ^A to % := L^(M'^) by 

L-.M^n-, i{y) = x{y;-) 

with Eq. 0. Then, it is easy to show that the embedding z. : ^ "H in fact satisfies condition 
Q of Proposition 2.1, where the almost complex structure Jy : Ty^A — TyA4 is given by 



Jy(q,p,A,B,<j},5^ 

= (B-\Aq-p), {AB-^A + B)q-AB'^p, -B, A, p'^B-~\Aq - p) - 6, -p-q + f), 
and hence M is symplectic. 
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Note that the variable 6 is essential in the symplectic formulation. We have 

,2 .li^ 



M{B,6) := \\x{y;-W 



detB 



exp 



25 



(9) 



and so we may eliminate 5 by solving ||x|| = 1 for 5 and substituting it back into Eq. ([7]) to normalize 
it. However, without 6, the manifold M is odd- dimensional and hence cannot he symplectic. More 
specifically, the variable 6 plays the role of incorporating the phase variable 4> into the symplectic 
setting. 

Remark 3.1. As we shall see later, Af{B, 6) = ||x||^ is essentially the conserved quantity (momentum 
map) corresponding to a symmetry in the system (by Noether's theorem). Normalization is intro- 
duced as the restriction of x to the level set ||x|| = 1 of the conserved quantity, i.e., x is normalized 
on the invariant submanifold of A4 defined by ||x|| = 1- Furthermore, this setup naturally fits into 
the setting of symplectic reduction and reconstruction as we shall see in Sections |4] and [6j 

3.2. Symplectic Gaussian Wave Packet Dynamics. We may now calculate the symplectic 
one- form Qm: Eq. Q, explicitly as 



:= L*Q = N{B, 5) ( Pi dg* - - iT{B-^dA) - d0 



(10) 



and hence also the symplectic form on M: 
Q,M ■= — d6^ 



: Af{B, 6) { dq' A dpi - ^dq' A tiiB'^dB) - ^dq' A dS 



2p^ 



h 



+ o (2^rfe'^H' + s^'sr,')dA, A dB, 



M 



1 2 1 

+ - [tT{B-UA) A d(5 - tr{B-^dB) A d(P] + -d(P A d6 \ . (11) 



On the other hand, the Hamiltonian becomes 

^2 



(12) 



where {V) {q,B, 6) is the expectation value of the potential V for the above wave function x^ i-C- 



(V) {q,B,6) :=exp(^-^) F(x) exp 

and {V){q, B) is a normalized version of it: 

{V)iq,B,5) ' 



(x — q)^B{x — q) 



dx 



{V){q,B) :- 



' detB 



M{B,6) 

In what follows, for any function A{x), we write 

(A) 



(A) :- 



M{B,5) 



V{x) exp 



(x - qfB{x 



dx. 



(13) 



\X\ 



\x\\ 



Note that if x is normalized, i.e., N{B,5) = \\x\\^ = Ij then {A) = {A); in particular {V) = {V). 
Now, the main result in this section is the following: 
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Theorem 3.2. The Hamiltonian system ixn^M = d-ff with the above symplectic form (11) and 
Hamiltonian (12) gives Heller's equations 1161 \17^ (see also Lubich \24\ Section II. 4-1]): 

(14) 



^, p = _(^s7V), A = --{A^ -13^)-{V^V), B = --{AB + BA), 
mm m 



{V)-^trB + ^ tr(B-\V^V)) , S = ^trA, 



2m 2m 4 V / 2m 

where V'^V is the d x d matrix defined by 

Proof. Calculation of ixn^M is straightforward, whereas that of dH is somewhat tedious: Note 
first that the derivatives of the potential term {V){q,B) are rewritten as foUows using integration 
by parts: 

As a result, we have 

dH = Miq, B) ( (W) • dg + ^ • dp + tr \(AB-'^ + B'^A) dA] 
\ m Am 

h 

+ -tr 
4 



-{Id - B-^A^B-^) - I HB-^ - B-'^{V^V)B- 
m n 



dfi^ -^Hd6 



where Id is the identity matrix of size d and H is what later appears as the reduced Hamiltonian 
in Eq. ([iS]): 

H:=^ + ^tr[B-\A^ + B^)] +W)iQ,B). □ 
2m 4m 

Remark 3.3. The original formulation of Heller |16j (see also Lee and Heller [21]) is not from a 
Hamiltonian/symplectic point of view and does not involves expectation values {V) etc. The above 
Hamiltonian formulation with the expectation values seems to be originally due to Faou and Lubich 
dU] and Lubich ^ Section IL4]. 

Remark 3.4. Writing C = A + iB, the above equations for A and B are combined into the following 
single equation: 

C = --C^ -lyW). (15) 
m 

4. Momentum Map, Normalization, and Symplectic Reduction 



The previous section showed that the symplectic structure for the semiclassical dynamics ( 14 ) is 
inherited from the one for the Schrodinger equation as its pull-back by the inclusion l : A4 ^ H. 
In this section, we show that the semiclassical dynamics also inherits the phase symmetry and 
the corresponding momentum map from the (full) quantum dynamics, and thus we may perform 
symplectic reduction, as is done for the Schrodinger equation in Marsden et al. |28i Section 5A] 
and Marsden [251 Section 6.3]. 

4.1. Geometry of Quantum Mechanics. Consider the S^-action : xH ^ H on the Hilbert 
space n = L2(M^) defined by 

The corresponding momentum map J : H ^ so(2)* = M, where we identified with S0{2), is 
given by (see, e.g., Marsden [25^ Section 6.3]) 

JW = -h\ 
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The expectation value of the Hamiltonian (H) is invariant under this action, and hence Noether's 
theorem imphes that the norm ||'0|| is conserved along the flow of the Schrodinger equation. In 
particular, the level set at the value —h gives the unit sphere S{T-l) in the Hilbert space T-L, i.e., the 
set of normalized wave functions: 

J-^{-H):={^€n\ 11^11= 1}=:S(?^). 

Since S-*^ is Abelian, the projective Hilbert space F{T-l) = J^^(— ^)/S^ = S{'H)/S^ is the reduced 
space in Marsden-Weinstein reduction [27] and hence is symplectic, i.e., with the inclusion i/j and 
projection vr^j defined by 



we have the symplectic form 0, on ¥(71) such that 

We may then reduce the dynamics to P(^). Note that the geometric phase (Aharonov-Anandan 
phase [1]) arises naturally as a reconstruction phase, as shown in Marsden et al. [28i Section 5A] 
and Marsden Section 6.3]. 

4.2. Geometry of Gaussian Wave Packet Dynamics. The geometry and dynamics in A4 
inherit this setting as follows: Define an S^-action <I> : x — t- on the manifold A4 by 

<^>e:M ^ M; {q,p,A,B,(j),5) ^ {a,p,A,B,(p + h9,6). 

Then, it is clear that the diagram below commutes, and hence <I> is the S^-action on A4 induced by 
the action ^' on 



M 



M 



The infinitesimal generator of the action with ^ £ so (2) 

d 



IS 



hi 



d 



e=0 



is defined by the condition 



The corresponding momentum map J_a4 : Ai — t- so(2)* S 

{^M{y),i) = {QM{y),iM{y)) 

for any ^ G so (2) and hence 

iM{y) = -^^{B,5). 

Thus, we see that J_a4 = J o t or Jx(y) = J(x(?/))- 

Now, the Hamiltonian : — )• M is invariant under the action, and hence again by Noether's 
theorem, J_a/( is conserved along the flow of Xh-, i.e., a level set of J_a/( is an invariant submanifold 
of the dynamics Xh- In particular, on the level set 

^m{-K) := {y e M \ iM{y) = -fi] , 

we have Af{B, 6) = 1 and thus, by Eq. ([9]), the Gaussian wave packet function x is normalized, i.e. 
11x11 = 1, and we may write 

'detB^'/' 



exp 



-{x - qy {A + iB){x - q) +p- {x - q) + 
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by eliminating the variable S as alluded in Section 3.1 Ignoring the phase factor e'''^^^ in the above 
expression corresponds to taking the equivalence class defined by the S^-action, and so the wave 
function 



detB 



[TTh) 



1/4 



exp 



1 



X — q)^{A + iB){x — q) + p ■ (x 



may be thought of as a representative for the equivalence class [xlj-i(_rj)] ™ the projective Hilbert 
space Fin). ^ 



Theorem 4.1 (Reduction of Gaussian wave packet dynamics). Heller's equations (14) on Ai is 
reduced by the above -symmetry to the Hamiltonian system 

ix^Tin = dH (16) 

defined on 



with the reduced symplectic form 
and the reduced Hamiltonian 



dq' A dpi + -B7^'B^/dAj A dBki 



H 



As a result, Eq. (16) gives the reduced set of Heller's equations: 

P . ; 1 



m 



p 



A = {A' 

m 



B' 



(V2y), 



B 



(17) 
(18) 

~{AB + BA). (19) 
m 



A few remarks are in order before the proof: 



term proportional to h. The quantum term is in fact essentially the imaginary part of the Hermitian 
metric 



Remark 4.2. Note that the reduced symplectic form Qfi is much simpler than the original one 
in Eq. (11); it consists of the canonical symplectic form of classical mechanics and a "quantum" 

itian 

(20) 
(21) 



5Sd 



ti(B-^dC B-^dC 



BT^}Bl^:^dCki ® dC^j 



on the Siegel upper plane [36j, i.e., 



lmgj,, = -BT^'Bi^dA, ^dBuu 



and this gives a symplectic structure on the Siegel upper plane 



Remark 4.3. Again, we may replace the last two equations of (19) by the succinct form 

1 



C 



m 



with C = A + iB. 



Proof of Theorem 4-1 > A simple application of Marsden-Weinstein reduction \27\ (see also Marsden 
et al. |29l Sections 1.1 and 1.2]). In fact, all the geometric ingredients necessary for the reduction 
are inherited from the (full) quantum dynamics as follows: Define the inclusion 



in : Jj{-h) ^M, 



the quotient map 



and also another inclusion 



[L]:Mn^n'H); [y] ^ [xiv)], 



10 
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where [ • ] stands for the equivalence classes defined by the S^-actions ^' and Then, the diagram 
below commutes and shows how the geometric structures are pulled back to the semiclassical side. 



M ^ 



Figure [T] gives a schematic of the inheritance. 




Figure 1 . Geometry of Gaussian wave packet dynamics: The geometric structures 
necessary for symplectic reduction of semiclassical dynamics on M are inherited 
from the full quantum dynamics in % as pull-backs by inclusions. 



The level set 3j^{—h) is defined by M{B,6) = 1, and so one may eliminate 6 (see Eq. ([9])) to 
write 

J^H-h) = X X §1 = {{q,p,A,B,^)}, 

and therefore the Marsden-Weinstein quotient is given by 



Then, the reduced symplectic form (17) follows from coordinate calculations using its defining 
relation 



We also have the reduced Hamiltonian H : Aiji — )• M, which appeared earlier in Eq. ( 18 ), uniquely 
defined by 

due to the S^-invariance of the original Hamiltonian H. 

Then, the Hamiltonian dynamics ixh^M = diJ on A4 is reduced to the Hamiltonian dynamics 
ix— rifj = dH on the reduced space Airi- D 
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5. Spherical Gaussian Wave Packet Dynamics 

This section is a brief detour into a simple special case of Gaussian wave packet dynamics that 
assumes that the wave packet is "spherical", i.e., A = aid and B = bid with Id being the identity 
matrix of size d; hence we replace the Siegel upper plane by Si even if d / 1. We also introduce 
the Darboux coordinates for the resulting semiclassical dynamics; they will be later exploited in 
the harmonic oscillator example in Section [8] to find the action-angle coordinates. 

5.1. Spherical Gaussian Wave Packet Dynamics. Setting A = aid and B = bid in Eq. ([T]) 

gives the "spherical" Gaussian wave packet, i.e., 

^(a + ib)\x - q\'^ + p ■ {x - q) + {(f) + i5) 



xiy;x 

The manifold M is now 



exp 



M 



T*R'^ X Si X 



X S\ 



Note that the Siegel upper plane Si = {a + i6 G C | 6 > 0} is literally the upper plane of C. The 
manifold Ai is {2d + 4)-dimensional, and is parametrized by 

y ■= {q,p,a,b,4>,S). 



The symplectic one- form 0x, Eq. (10), now becomes 



:= i*e = Af{b, 6) ( Pi dq' - ^ da - d</. 



with 



TTh\ 



d/2 



M{b,6) : 

and hence the symplectic form on Ai is 
0_A^ = J\f{b, 6) dq^ A dp 

d{d + 2)h 



exp 



25 



^dq' A d6 - ^dq' A d6 
2b n 



+ 



da A d6 + 7^(do A d(5 - d6 A d(/>) + ^d(/> A d5 



852 ' 2b 

which is given by Faou and Lubich [TU] (see also Lubich [2T, Section II.4]). 
On the other hand, the Hamiltonian : — )■ M, Eq. ([s]), is given by 



H = M{b,6) 
= M{b,6) 



2m 

1 
2m 



2 rr-a^ + b 

p' + dn — 
p +dh^^ 



+ {V) {q,b,6) 



(22) 



where 



and 



(l') {q,b,S) -.= {v,Vx> = 



exp 



25 



V{x) exp ( 



6, 



xd/2 r. 



M{b,6) 

Hence, as shown in |10j . the Hamiltonian system ([s]), i.e., 
with 



V{x) exp( ~^\x 



dx, 



dx. 



(23) 



d 



d 



d ; d 



oq^ opi oa ob 



— 6— 
d(h 85 
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gives the spherical version of the equations of HeUer [16] : 



^, p = -vy, a = -A^^)^ b = ' 

m mam |-2^-j 

; p"^ TTTT dh , h -r-r—- ■ dh 

<p = - (V) b+—(AV), 6= — a. 

^ 2m ^ ' 2m 46^ " 2m 



We may apply the symplectic reduction in Theorem |4.1| to obtain the following reduced symplectic 
form on Adu: 

= dq' A dpi + -r^da A db. 
4o^ 



The reduced Hamiltonian ( 18 ) is now 



H=^ + dh-±- + {V){q,b), 
2m 4m b 



and the reduced equations (19) become 

q=^, p=-{VV), a = -Ay, b= . (25) 

m mam 

5.2. Darboux Coordinates. Let us define the new coordinate system 

{q\r,^,Pi,Pr,p^) ■.= (^q\ -c^, Af{b, 5) p, ^a,Af{b,5)j . (26) 

Then, the symplectic form takes the canonical form 

flM = d(7* A dpi + dr A dpr + dip A dp^, 

and thus the above coordinates are the Darboux coordinates. Hence, the Hamiltonian system (25) 
is transformed to the following canonical form: 

,■ dH dH dH 

r = 7^, f 



dpi dpr dp^ 

■ _ dH . _ dH . _ dH 

Pi — o~T' P^ — o ' Pv — o ■ 
oq^ or 0(p 



The Darboux coordinates ( |26[ ) for Ad induce those for i-e-) noting p = p, we have the 

Darboux coordinates {q^,Pi,r,pr) for Aiji with 

{r,pr) = —i-,a 



that is, the reduced symplectic form ^Ifi takes the canonical form 

Ufi = dq' Adpi + dr Adpr. (27) 



Therefore, the reduced dynamics (16) with 



^ dq^ ^ dpi dr ^ dpr 
is written as canonical Hamilton's equations: 

, dH dH dH dH 

opi dpr dq^ dr 

Remark 5.1. That the variables b~^ and a are essentially canonically conjugate was pointed out by 
Littlejohn |23| and Simon et al. [38]. See also Broeckhove et al. ^TJ and Pattanayak and Schieve 
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6. Reconstruction — Dynamic and Geometric Phases 



6.1. Theory of Reconstruction. As described in Section 4.2 , the Gaussian wave packet dynamics 
Xh defined by Q in M may be reduced to the Hamihonian dynamics Xj^ defined by ( |16[ ) in the 
reduced symplectic manifold Airi '■= -^m{~^)/^^- Now let c{t) be an integral curve of the reduced 
dynamics Xjj, i.e., c{t) = Xjj{c{t)). Then, the curve c{t) is the projection of an integral curve 
c{t) of the full dynamics Xh on J^(— i.e., vr^i o c{t) = c{t). Then, a natural question to ask is: 
Given the reduced dynamics c{t), is it possible to construct the full dynamics c{t)? The theory of 
reconstruction (Marsden et al. [28]; see also Marsden [25j Chapter 6]) provides an answer to the 
question, and the so-called dynamic and geometric phases arise naturally when reconstructing the 
full dynamics from the geometric point of view. 

6.2. Dynamic Phase. For the full quantum dynamics with the Schrodinger equation, we may 
define a principal connection form £/ : TJ~^(—h) — )• so(2) on the principal bundle J~^{—h) — )• F(T-l) 
as follows (see Simon |37] and Montgomery [33j Section 13.1]): 

= Im(V',dV') |r^j-i(_ro; {£/ {ip) , v^) = Im {4) , v^) for e T^J'^{-h), 

that is, £/ is —Q/h restricted to J~^{—h). Since HV'II^ = = 1 for V S J~^(— ^), we have 

{dip,ip) + {ip.d-ip) = 0, and thus 

jz/ijP) = -i {^,d^) . 

This induces the principal connection form ■ ~^ so(2) on the principal bundle 

Jj^{—h) —7- Aifi that is given as follows: 

which, for the spherical case, reduces to 

.<2/m = - T Pi dg* + -J^ da = ^ (-d^p - pi dq' + r dpr) . (29) 

Now, let yo be a point in JjJ^{—h) and d{t) be the horizontal lift of the curve c{t) such that 
d(0) = yo, i.e., the curve defined uniquely by tth o d{t) = c{t) and d{0) = yo with 

.s^M{d{t)) ■ d{t) = 0. (30) 

Then, since the full dynamics c{t) satisfies vr/j o c(t) = c(t), we have tt^^ o c{t) = tt/j o d{t), and thus 
there exists a curve g{t) in such that c{t) = g{t)d{t). By the Reconstruction Theorem ( |28| 
Section 2A] and j25l Section 6.2]), the curve g{t) in is given by 

g{t) = exp(i [ as)ds) , (31) 



'M := t*£/ = i:dcP - ^pidq' + ^tv{B-'dA), (28) 



where 

m ■■= ^M{d{t)) ■ XH{d{t)) 

is a curve in so (2) = M. It is straightforward to see, from Eqs. (12), (25), and (29), that 



H{d{t)) H{c{t)) E 

^^'^ = — tr- = — tr- = -T ^^^^ 

where the the second equality follows from the S^-invariance of the Hamiltonian H; the last equality 
follows since c{t) is an integral curve of Xh, and so the Hamiltonian H is constant along c(t), and 
its value is determined by the initial condition E := H{c{0)). Therefore, we obtain 

g{t) = exp(-'-Et 
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which is compatible with the result for the full quantum dynamics (see, e.g., Montgomery |33| 
Section 13.2]). 

In particular, if the curve c{t) is closed with period T, i.e., c(0) = c(T), then the dynamic phase 
^dyn G achieved over the period T is given by 

gdyn = exp^^A(/)dyn^ = exp^-^E'T^ . 
where A(/>dyn is the change in the angle variable (f> in the coordinates for M (see also Eq. ([T])): 

A0dyn = -^T. 



6.3. Geometric Phase. The curvature of the principal connection form (28) is given by 



= d£/M = \ {^q' A dpi + ^B-^^Blj^dAij A dBki^ , 
and, for the spherical case, we have 

= ^ (^dg* A dpi + A db^ . 

Therefore, its reduced curvature form, i.e., viewed as a two- form on Ain, becomes 

= l^h- 
n 

Suppose that the curve of the reduced dynamics on A^n is closed with period T, i.e., c(0) = c(T) 
for some T > 0. Then, the geometric phase (holonomy) 5geom G achieved over the period T is 
defined by 

d{T) = ggeomd{0). 

Let D be any two-dimensional submanifold of Airi whose boundary is the curve c([0,T)); then 
the geometric phase is given by the following reconstruction phase (see, e.g., Marsden et al. |28| 
Corollary 4.2]): 

5geom = exp^^A(/)geom^ = exp^-i j j ^ti^ = exp^-^ jj G S\ 

where A(/>geom is the change in the angle variable (j): 

A(/)geom = - jj ^h- (33) 

This generalizes the result of Anandan [31 HI [5] , which was derived for the frozen Gaussian wave 
packet, i.e., the spherical case with a and b being constant. 

Notice that we derived the above formula as a reconstruction of the Hamiltonian dynamics in 
M; namely, we have incorporated the phase variable (j) (accompanied by 5) into the expression of 
the Gaussian wave packet ([T]) to write the full dynamics in as a Hamiltonian system (see the 
discussion just above Remark |3.1[ ), and the reconstruction of the dynamics on Jv[ from the reduced 
dynamics on A^/j gave rise to the geometric phase. This gives a natural geometric account (and 
generalization) of the somewhat ad- hoc calculations performed in [SHS] . 
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6.4. Total Phase. Combining the dynamic and geometric phases, we obtain the total phase change 
over the period T: 



9total 



exp( -A(/>total 



5dyn ■ 5geom — exp 



ET + 



n, 



D 



or 



A(;^total = Ac/fdyn + A(/), 



geom 



-ET 



D 



which is similar to the rigid body phase of Montgomery [32] (see also Hannay |15j . Anandan [5], 
and Levi [22j). Noting that the phase factor in ([7]) is e*"^/^', it is convenient to rewrite the result as 

1 

h' 



A 



'total \ 



-ET 



D 



(34) 



Note that we made an assumption that the reduced dynamics on Mfi defined by Xjj is periodic 
with period T. In Section 



8.3 



below, we will show that such a periodic orbit in in fact exists 
for the semiclassical harmonic oscillator and calculate the explicit expression for the total phase. 

If the reduced dynamics is not periodic, we do not have a simple formula for the phase change as 
above. However, one may still obtain an expression for the phase factor (f> in terms of the reduced 
solution c{t) = {q{t),p{t),A{t),B{t)) defined by Eq. Let us write 

dit) = {qit),pit),Ait),B{t),m,m), 
Since c{t) = g{t) d{t) with g{t) given by Eq. ([31 



c{t) = {q{t),p{t),A{t),B{t),(t>{t),6{t)). 



and so, using the expression for ^{t) in Eq. (32), 



-H{d{t)) + ^{t). 



Now, the horizontal lift equation ([30|) gives 



m 



+ ^ti\B-\A^-B^)\+^t4B-\VW) 
Am -'4 V 



where we used the reduced equations (19). As a result, by using the expression for the Hamilton- 
ian (12) and noting that N{B, 6) = 1 here, we obtain, 

h h 



— - {V) - — tj:B+'-^ tr( B-\V^V) 



which is simply the equation for (p in the full dynamics (14). 



7. Asymptotic Evaluation of the Potential 



One obstacle in practical applications of Heller's equations, Eq. (14) or (19), is the evaluation 
of the potential terms (V), (VV), and (AT/), which are generally given by complicated integrals 
(see Eqs. (13) and ([23|)). If the potential V{x) is given as a simple polynomial, one may reduce the 
integrals to Gaussian integrals and obtain closed forms of them exactly; this is particularly easy for 
the spherical case (see Eq. ([23[)). However, one rarely has such a simple potential V{x) in problems 
of interest in chemical physics, and thus there is a need to approximate the potential terms. 
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As mentioned in Remark 3.3 , Heller's formulation does not involve these averaged potential terms, 
but from our perspective, it can be interpreted as adopting the following simple approximations of 
the expectation values: 



{V)iq, B) ^ V{q), {VV)iq, B) ^ Wig), {AV)iq, B) ^ AViq). 

Notice, however, that this approximation neglects non-classical effects coming from B altogether, 
and seems to be too crude for a semiclassical model. 

Instead, we apply Laplace's method to the integral in the potential term {V) to obtain an 
asymptotic expansion of it. As we shall see later, this also results in an asymptotic expansion of 
the Hamiltonian H, Eq. (12), and then our Hamiltonian/symplectic viewpoint provides a correction 
term to the formulation by Heller |16j and Lee and Heller j21j . 



7.1. Non-spherical Case. The key observation here is that the potential term {V), Eq. (13) or 
(23), is given as a typical integral to which one applies Laplace's method for asymptotic evaluation 
of integrals, i.e., we have 



where 

Fn{q,B):= I e^(^)/'^y(a;)(ix (35) 



with 

R{x) = -[x-qYB{x-q). 
Now, an asymptotic evaluation of the integral F}i{q, B) gives us the following: 

Proposition 7.1. // the potential V{x) is a smooth function such that e'^^'^^^^^V{x) is square 
integrable in M'^ for some a G [0, 1), then the potential term {V) has the asymptotic expansion 

oo 

'{V){q,B) ^^Cn{q,B)h'' as h^O, (36a) 

n=0 

where 

Pk all even 

and 

g,{0 := D^V{Qi) = g^Pi , , ggg. ^(20, (36c) 

with V{^) := V{q + bi, . . . b^ are the eigenvalues of B, and Q is the orthogonal matrix such that 

BQ=Qdiag{bi,...,bd), 
i.e., each of its columns is an eigenvector ofB. 

Proof. The asymptotic expansion follows from a standard result of Laplace's method (see, e.g.. 
Miller [311 Section 3.7]) applied to the integral Ffi{q,B) in Eq. (35) restricted to a neighborhood of 
the point x = q. Hence, we need an estimate of the contribution from the remaining part of the 
integral to justify the expansion. See Appendix \K\ for this estimate. □ 

In particular, we can rewrite the first two terms more explicitly: 

ly){q, B) = V{q) + ^ tr [B~^V^V{q)\ + 0{h^) as h ^ Q. (37) 
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Therefore, the Hamiltonian (12) becomes, as h 

^2 



H = M{B, + V{q) + -P- tilB-HA^ + B^)] + ^ tr(B~'V^V{q)) + OiH^ 

I Zm Am ^ -'4 

We may then neglect the second order term 0{h?) to obtain an approximate Hamiltonian 



~ i/i := M{B, 5) 



2m 4 



B- 



A^ + B^ 



m 



Then, the Hamiltonian system ixjj_^^M = dHi gives the following approximation to Eq. (14): 



P 

q = 

m 

A = --{A 

m 



P 



d_ 
dq 

-B'^)-V^V{q) 

— -V( )- — 
2m 2m 



V{q) + -tT{B~'V^V{q)) 



B = 
tiB, 6 



~{AB + BA), 
m 



(38) 



2m 



tiA. 



Notice a slight difference from those equations obtained in Heller [16] and Lee and Heller |21j : The 
second equation above has a semiclassical correction proportional to h, whereas those in 116\. \21^ are 
missing this term. Furthermore, since the correction term generally depends on B, the equations 
for q and p are not decoupled as in Heller fl^ . Therefore, it is crucial to formulate the whole 
system — as opposed to those for q and p only — as a Hamiltonian system. 



Remark 7.2. If the potential V{x) is quadratic, then the asymptotic expansion (36) terminates at 
the second term, i.e., c„ = for n >2, and becomes exact. Hence H = Hi and so Eqs. (14) and 
(38) are equivalent. Moreover, since V'^V{q) is now constant, the second equation in (38) reduces 



to the canonical one, and hence Eq. (38) reduces to those of Heller [16j and Lee and Heller 

The reduced system ix = dHi with the 



One may reduce Eq. (38) just as in Theorem 
reduced approximate Hamiltonian 

„2 



4.1 



Hi :-- 



^ + V{q) + ^tr 
2m 4 



B- 



A^ + B^ 



m 



+ W(g) 



gives the first four equations of ( 38 ) . Notice that the Hamiltonian is split into the classical one and 
a semiclassical correction proportional to h. 



7.2. Spherical Case. The asymptotic expansion for the spherical model from Section [5] follows 
easily from Eq. (37): Setting B = hid gives 



{V){q,B) = V{q) + -^V{q) + 0{h^) 



as 



Note that higher-order terms are easy to calculate for the spherical case, because B = hid implies 
that hk = h for k = 1, . . . , d and Q = Id- Now, Eq. (38) becomes 
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8. Example: Semiclassical Harmonic Oscillator 

In this section, we illustrate the theory developed so far using a simple one-dimensional harmonic 
oscillator. For this special case, the system (25 ) is easily integrable as shown by Heller [16]; however, 
we approach the problem from a more Hamiltonian perspective. Namely, we first find the action- 
angle coordinates for the reduced system using the Darboux coordinates from Section 5.2 As we 
shall see later, the action-angle coordinates give an insight into the periodic motion of the system, 
and facilitates our calculation of the geometric phase. 

8.1. The Hamilton Jacobi Equation and Separation of Variables. Consider the one-dimensional 
harmonic oscillator, i.e., d = 1 and 

1 



V{x) 



2 2 

-moj X . 



Note that for the one-dimensional case, the non-spherical wave packet reduces to the spherical one. 
Then, the potential term is easily calculated to giv^ 



V{q,b) 



and so the Hamiltonian ( 22 ) is 



1 
2m 



+ dh- 



46 

2 , 52 



26 



+ 



muj 



+ 



26 



or, using the Darboux coordinates defined in Eq. (26) 

o2 



H 



2m p^ 



muj"^ 2 2 2 mu)"^ 
2 ^ m 2 8mr 



Then, the reduced Hamiltonian ( 18 ) becomes 

a2 + 62 



H 



1 

2m 
1 



p'^ + dh 



26 



+ 



moj 



9 h 



— {p'' + 4.rpl)+—{q^ + r) + 



8mr ' 



which also follows from Eq. (39) with p^ = 1. 

The Hamilton-Jacobi equation for the reduced dynamics 

dW dW 
dq ' dr 

with the ansatz W{q,r) = Wg{q) + Wr{r) gives 



H[q,r, 



E 



1 fdW, 
2m \ dq 



+ 



2r fdWr 



muj , n 
+ ((?'+r) + 



m \ dr 

Hence, by separation of variables, we obtain 

1 fdW, 
2m \ dq 

where Ei and Er are constants such that Ei + Er = E. Thus, 



8mr 



2 2 
moj n 



^1 



2r (dWr 
dr 



m 



mu! 

_l ^ . 

2 



E. 



8mr 



Er 



dWg 

dq 



±y/2mEi — m^uj^cp'^ 



dWr 

dr 



mui 



a 
r 



2r2' 



(39) 



Since V{x) is quadratic, the asymptotic expansion (36 1 is exact, i.e., c„ = for n > 2, and so (36 1 gives the same 



result. 
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where 

a := ^, L 



and we assume that L < a/\/2 is satisfied. 



8.2. Action— Angle Coordinates. The above solution of the Hamilton-Jacobi equation gives 
rise to the canonical coordinate transformation to the action-angle coordinates, i.e., {q,r,p,pr) i— t- 
(6*1, 9r, h,Ir)- 

The first pair of action-angle coordinates {9i,Ii) are those for the classical harmonic oscillator: 
Let 7i be the curve (clockwise orientation) on the q-p plane defined by p'^ ={dWg/dq) , i.e., 

1 2 muP' o „ 
2m 2 

which is an ellipse whose semi-major and semi-minor axes are y^2Ei/m/(jj and \j2'mE\. Therefore, 
the first action variable I\ is given by Stokes' theorem as follows: 

h = — f pdq = — dpAdq= — = - —p H — 

2tt J^^ 27r uj L0\2m 2 

where Ai is the area inside the ellipse (with the orientation compatible with that of 71), i.e., 
dAi = 71; hence the surface integral is the area of the ellipse. The angle variable 61 is then 

Interestingly, the second pair of action-angle coordinates {9r,Ir) is essentially the same as those 
for the radial part of the planar Kepler problem (see, e.g., Jose and Saletan [20^ Example 6.4 on 
p. 318]). Let 7r be the curve (clockwise orientation) on the r-pr plane defined by p"^ ={dWr/dr) , 
i.e., 

2r n moj'^ It? 



[ dWq f r- , -if'mujq\ 

/ dq= \ y 2m uoli — m'^uj'^q'^ dq = tan I 

J dq J \ P J 



— pl+'^r+ — =Er, 
m 2 6mr 



or 



moj a 

^^ = ^^V-' + 7-2^- 

Setting pr = yields r = r± := (a it \/ o? — 2L'^)/2. Then, the action variable Ir is calculated as 
follows: 



Er h r 2 "^'^ ^'^ ^ 

= — — — = p -\- r + — — . 

2uj 4 muj ^ 4 16m ur 4 

The angle variable Or is then given by 



d fdWr ^ ^ _1 

dr = tan 



dir J dr 



4r'^{m'^u}'^ - 4p2) - fi^ 
16m UJ r^pr 



The (reduced) Hamiltonian H is then written in terms of the action variables as follows: 

H = [ll+ 2Ir + ^ j UJ. 

Then, the reduced dynamics on Aifi is written as 
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and Ii and Ir are constant. Therefore, the reduced dynamics is now transformed to a periodic flow 
on the torus = x = {((9i, 6',.)}. 

8.3. Calculation of Geometric Phase. Recah, from Section [6| that we may calculate the geo- 
metric phase achieved by a periodic motion of the reduced dynamics on A4ri- The previous section 
revealed that the reduced dynamics is in fact periodic with period T = 2-k/u}; we have also obtained 
the curves traced by the periodic solution on the q-p and r-pr planes. These results enable us to 



calculate the geometric phase explicitly. First recall from (33) with (27) that we have 



A(/>geom = - {dq Adp + dr Adpr), 



ID 

where D is any two-dimensional submanifold in Aiji whose boundary is the periodic orbit T G Ain, 
i.e., the curve c : [0, T] — )• Aifi defined by the reduced dynamics. Then, the projections of the curve 
r to the q-p and r-pr planes are the curves 71 and j,. defined above, including the orientations 
(note that the clockwise orientations for 71 and 7,. coincide the direction of the dynamics on F). 
Therefore, we have 



A(/>geom = f pdq + 2 d) prdr, 

since the projection of F to the r-pr plane gives two cycles of for a single period T = 27r/uj. 
Using the expressions for p and pr from the above subsections, we obtain 

A0n;oom = irh = ET — irh, 

OJ 

which gives the following Aharonov-Anandan phase (note that the phase factor in ([T]) is e*'^/'*): 

. / 0gcom\ ET 



and hence the total phase change is given by, using Eq. (34), 

^ ^ '/'total ^ _ ^ ^ 4>djn + 0gcom ^ _ _^ 

This implies that the corresponding wave function (see Eq. ([7])) flips "upside down" (just like a 
falling cat!) after one period, i.e., 

ioy(r) = -ioy(O) or x{y{T),x) = -x{y{^),x). 



9. Geometry of the Hagedorn Wave Packet Dynamics 

Hagedorn [13j used a slightly different parametrization of the non-spherical Gaussian wave packet 
([7]) in his study of the asymptotic behavior of Gaussian wave packet dynamics in the semiclassical 
limit ^ — )■ 0. This formulation also leads to an elegant derivation by Hagedorn [H] of an orthonormal 
basis generated by the Gaussian wave packet with certain raising and lowering operators. Faou 
et al. [IT] exploited this orthonormal basis to develop a numerical method to solve the Schrodinger 
equation more efficiently than with the Fourier basis. 

In this section, we give a geometric account of Hagedorn's wave packet dynamics, building on 
the results by Lubich [24', Chapter V]. The main objective is to understand its relationship with 



Heller's equations ( 14 ) from the geometric point of view. 
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9.1. The Hagedorn Wave Packet Dynamics. The Hagedorn wave packet [13] (see also Lubich 
\24:\ Chapter V]) parametrizes elements C = A + iB in a different way; specifically, one writes 
C = PQ-^ in Eq. ^ to have 



where y := {q,p,Q,P,(l),6). Its norm is then 



x{y;x) = exp<{ - 



1. 
2^ 



+ p ■ {x - q) + {(j) + i6) 



M{Q,5) := My, ■)f = {nhr/^\detQ\exp(^-'^'^ , 
and JM{y) = —hM{Q, 6); hence the wave packet is normalized on the level set Jj^{—h) as follows: 



exp< 



(40) 



where we defined the new variable 



S 



arg(det Q) 



and the "ground state" of the Hagedorn wave packets ipo (Hagedorn [13] and Lubich 
ter V]): 



Chap- 



Mq,P,Q,P;x) :=(^/i)-'^/^detQ)-^/2gxp|^ _ _ ^) + p . _ | 



h 

Hagedorn |14j showed that Eq. ( |40| ) is an exact solution of the Schrodinger equation with a 
quadratic potential V{x) if the parameters {q,p, Q, P, S) satisfy 

P_ 
m ' 



Q 



p 

m ' 



q = 
p 



p=-VV{q), 



-VV(g)Q, 



(41) 



Note that the last equation says that the variable 5 is essentially the classical action integral, i.e.. 



S{t) = 3(0) + 



2m 



V{q{s)) ds 



In the rest of the section, we look into the geometry behind the dynamics of the Hagedorn wave 
packet. 

9.2. Geometry of the Siegel Upper Plane S^. As we saw above, Hagedorn's wave packet is 
a different parametrization of the d x d complex matrix C = ^ + i;B in Heller's wave packet Q. 
Recall that C belongs to the so-called Siegel upper plane defined as (see Eq. ([s])) 

Srf := ^A + iBe C^""^ \ A,Be Symrf(M), S > o} . 

The key to understanding the geometry of the Hagedorn wave packet dynamics is the fact that the 
Siegel upper plane E^^ is a homogeneous space. Specifically, we can show that (see Siegel |3Bj and 
also Folland [12 Section 4.5] and McDuff and Salamon [301 Exercise 2.28 on p. 48]) 

j:d = Spi2d,R)/U{d), 

where Sp{2d,M) is the symplectic group of degree 2d over real numbers and U{d) is the unitary 
group of degree d. In fact, consider the (left) action of Sp{2d,M) on defined by 

'A B 
C D 



^ : Sp{2d,R) X Sfl 



^ ^ ^{C + DZ){A + BZ)-\ 



(42) 
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This action is transitive: By choosing 

h - ^1 - Ol 1 _ r 

[C d\ - [A h\ [ B^/^\ ~ i3i/2 

which is easily shown to be symplectic, we have 

Now, the isotropy group of the element ild G is given by 

G M2d(M) I U^U + V^V = Id, U^V = V^U 



Sp{2d. 



U V 
-V u 



Sp{2d,R)nO{2d), 



where 0{2d) is the orthogonal group of degree 2d; however Sp{2d,'R.) fl 0{2d) is identified with 
U{d) as follows: 

/ : 5p(2(i,M)nO(2d) ^ C/(d); _^ ^ ^ U + iV. 

Hence 5p(2(i,R)i/^ ^ U{d) and thus ^ 5p(2d,M)/[/((i). Indeed, we may identify Sp{2d,m}, /U{d) 
with Srf by the following map: 

Sp{2d,M)/U{d) ^ Srf; [g]u(d) ^ ^^giUd), 

where [ ■]u{d) stands for a left coset of U{d) in Sp{2d, M); then this gives rise to the following explicit 
construction of the quotient map S'p(2d,M) — )■ Sp{2d,]R.)/U{d); 

TTu^d) ■■ Spi2d, R) ^ Sp{2d, R)/U{d) ^ Srf; ^g{ild). 

Therefore, we have the following diagram, which simply shows that the action ^ is indeed a left 
action: Note that the map Sp{2d, M) x Sp{2d, M) — >• Sp{2d, M) is the standard matrix multiplication. 



Sp{2d,] 



Sp(2d,M) 
> Srf 



5 



-> hg 



(43) 



As shown by Siegel [36] , the map ■ '^d ^ '^d is an isometry of the Hermitian metric ( 20 ) for any 
h S Sp{2d,M.) and therefore is symplectic with respect to the symplectic form (21). This suggests 



that the S^-component of the symplectic dynamics defined by Heller's equations ( 19 ) may be lifted 
to the symplectic group Sp{2d,M.). It turns out, as shown in the next subsection, that this is 
precisely what the Hagedorn parametrization does. 

9.3. The Hagedorn Parametrization. As pointed out by Lubich [24] Section V.l], the parametriza- 
tion of Tid by Hagedorn [13i, J4j fits into this geometric picture very well. Lubich |24[ Section V.l] 
shows that the symplectic group 5^(2^,1^) is written as follows: 

"ReQ ImQ" 
ReP ImP 

= {(Q,P) G Mrf(C) X MdiQ I Q^P - P^Q = 0, Q*P - P*Q = 2ild} , 

and also that Hagedorn's parametrization of is nothing but the explicit description of the 
quotient map 

[ReQ ImQ" 
ReP ImP 



Sp{2d, 



E M2d(M) I Q, P e Md{C), Q^P - P' Q = 0, Q*P - P*Q = 2ild 



vr 



U{d) ■ Sp{2d, 



or (Q,P) ^ PQ'\ 



(44) 



This map provides the connection between the dynamics of Heller and Hagedorn alluded above: 
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Proposition 9.1. The Tid- component of Heller's equations, i.e., 



C 



1 



m 



is the projection to the Siegel upper plane of a curve 

"ReQ(t) lmQ{t) 



(15) 



g{t) 



ReP(t) ImP(t) 



in the symplectic group Sp{2d, M) defined by 

9{t) = TeRgit){m) = m9{t) 

with ^{t) e 5p(2(i,M) being 



or equivalently. 



m ■■- 



Q 



Id/m 



P 

m' 



-{V^V)iq,Q)Q. 



Remark 9.2. The expectation value (V^V) originally depends on B, but since B = (QQ*) ^ (see, 



e.g., Lubich [Ml Lemma V.1.1 on p. 124]), it now instead depends on Q. If V{x) is quadratic, then 
V^y is constant and {V'^V){q,Q) = V'^V{q) and thus the above equations for Q and P reduce to 
those of Hagedorn, i.e., Eq. (41). 



Proof of Proposition 9. 1 . Let C{t) be a curve in defined by ( |15| ) . Since Sp{2d, M) acts transitively 
on Tid, there exists a corresponding curve h{t) in Sp{2d,R.) such that "^^^{0(0)) = C{t) and 
/i(0) = Id- Now, let qq £ Sp{2d,M.) be an element such that T^u{d){9o) = C(0), and define the curve 
g{t) := h(t)g(0). Then, clearly we have iTu(d) ° g{t) = C{t), i.e., the following diagram commutes as 
in (|43]). 

Sp{2d, M) go . > h{t)go 



Sp(2d. 



■^U{d) 



"hit) 



Let us then write 



m ■.= g{t)9{t)-^ = h{t)h{t)-\ 
which is in the Lie algebra sp(2(i,M); thus it takes the form 

with and .^21 both being symmetric; then it is easy to see that g = ^g gives 



Therefore, from (44) and the above expressions, we have 

Tg7:u(d){g) = PQ-^ - PQ-^QQ-^ 

= '^21 — CnC — C^ii — C.^i2C. 

where we also used the relation 7^u{d)°9{'t) = C{t), i.e., PQ~^ = C. However, since 7^u{d)°9{t) = C{t), 
we have T-Ku(^d){9) = which implies, using the above expression and Eq. (15), 



^21 - iiiC - C^ii - C^i2C 



m 



- (V2y). 
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Therefore, we may take 

1 



61 = 0, a2 = -/d, 6i = -(W). □ 

m 

Remark 9.3. One may recognize ( |15[ ) as an example of the matrix Riccati equation. Particularly, if 
(SI'^V){q,Q) is independent of Q (which is the case if V{x) is quadratic, as mentioned above) and 
q{t) is known, then it is a matrix Riccati equation. In fact, there is a similar geometric structure 
behind the matrix Riccati equation: As shown in Hermann and Martin [TH] and Doolin and Martin 
[S] , one considers the action of a general linear group on a Grassmannian using the linear fractional 



transformation of the form (42), and then performs virtually the same calculations as above to 



derive the matrix Riccati equation. 



We may now rewrite Heller's equations ( 14 ) using Hagedorn's parametrization: 



p p 



p = -{VV), Q=-, P = -{VW)Q. 
m m 

The dynamics is now on the manifold M = T*W^ x Sp{2d, M) x x M; its dimension is 2(f + 3d + 2, 
which is odd if d is odd. Thus, unfortunately, the above system of equations cannot be Hamiltonian 
on M in the symplectic sense for an arbitrary d. 



10. Conclusion 

We gave a symplectic-geometric account of Heller's semiclassical Gaussian wave packet dynamics 
that builds upon on a series of works by Lubich and his collaborators. Our point of view is helpful 
in understanding how semiclassical dynamics inherits the geometric structures of quantum dynam- 
ics. Particularly, the geometry behind the symplectic reduction and reconstruction of semiclassical 
dynamics is inherited from quantum dynamics in a natural way. We also exploited the geometry 
of the Siegel upper plane to interpret the Hagedorn parametrization of the Gaussian wave packet 
from the geometric point of view. Furthermore, we derived an asymptotic formula for the expected 
value of the potential to approximate the potential terms appearing in the system of equations for 
semiclassical dynamics. The asymptotic formula not only naturally generalizes Heller's approxi- 
mation but also indicates that it is crucial to couple the equations for the classical position and 
momentum variables q and p with those of the other quantum variables, thereby justifying our 
point of view of regarding the whole system as a Hamiltonian system. 
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Appendix A. Proof of Proposition 17.11 
First take the ball Bg{q) and then split the integral Ff-^ in Eq. (|35|) as follows: 



g"v-y/" v{x) dx + 



,Ri^)/f^Vix)dx. 



Then, since V is assumed to be smooth, we obtain the asymptotic expansion (36) by applying the 
standard result of Laplace's method (see, e.g.. Miller [311 Section 3.7]) to the first term. We need 
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an estimate of the second term to justify the expansion. Introducing the variable ^ := x — q and 
defining 

R{^):=R{q + ^) = -eB^ and V{i) ■.= V{q + i), 

we have 



where we used the Cauchy-Schwarz inequahty. The second term is bounded by assumption. To 
evaluate the first term, we introduce the new variable rj = Q^^; then the exponent simplifies to 

d 

2(1 - u)R{Qn) = -2(1 - a) rj^O'BQri = - ^ hvl 

k=l 

where we set Pk '■= 2(1 — cr)6fe, which is positive. Thus, we have 



d 



where C^^^ is the hypercube defined by 

Ce/Vd ■= \Vk\ < ^ for fe = l,...,c?|, 

which is clearly contained in ^£(0). Writing ea ■= s/s/d for shorthand, the Cauchy-Schwarz 
inequality gives 



roo 



/ ft \ 3/4 



Therefore, 



8ey \2{1 -a)bj 

d/i / f^d \3/4 . 2/ 
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as h for any real p, since the above exponential term is dominated by any real power of h. 
Therefore, 



/ 

Jm. 



e^^^y^V{x)dx = o{tf) as ^^0 

<i\BM 



for any real p as well, and so the above integral has no contribution to the asymptotic expansion. 
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